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ABSTRACT ^“536 

\ 

A visco-plastic flow theory for the solution of hyper- 
velocity perforation of thin plates has been investigated. A 
short time analytical solution including the effects of target 
material yield strength on perforation radii and ballistic 
limits was obtained. Expressions were derived for the veloc- 
ities, displacements, strains, and strain rates present in the 
target material during hypervelocity impact. Pertinent param- 
eters were varied and a comparison was made with the solution 
in which the yield strength was taken to be zero. The effect 
of target material yield strength was not found to be signif- 
icant and resulted in little variation in the computed perfora- 
tion radii. 

Ballistic limit studies were also performed in which 
minimum plate thicknesses necessary to prevent perforation were 
established for projectiles of various masses, radii, and 
impacting velocities. Yield strength considerations were 
found to be of secondary importance; the more critical effect 
being produced by the projectile aspect ratio at impact. 
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INTRODUCTION 


An analytical solution to the problem of a rigid circular 
cylindrical projectile of finite mass impacting upon an 
infinite plate is presented in which yield strength considera- 
tions have been included in the governing equations. The 
results of reference 1 have been extended and, in addition, 
ballistic l imi ts of thin plates have been determined. For the 
sake of completeness, the solution which was presented in ref- 
erence 1 has been repeated in the appendix. 

In references 2 and 3 an analytical solution was obtained 
for a problem similar to that of reference 1 except that the 
yield strength of the target plate was taken to be zero in the 
governing equations. A nonzero yield strength of the target 
material was introduced subsequently only to establish a sep- 
aration criterion for perforation. An approximate two- 
dimensional analysis of the plate perforation problem (as 
stated in ref. 2 ) is presented in reference 4. This solution 
indicates that the radial velocities are about one order of 
magnitude less than the axial velocities and thus helps to 
justify the use of a one-dimensional approach. 

An extension of the work presented in references 2 and 3 
is given in reference 5, in which an attempt was made to 
include the effects of target material yield strength in the 
governing differential equation and associated boundary con- 
ditions. The results of the analysis are given in finite 
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series form and a parametric study is presented in which the 
effects of variations in the pertinent parameters on the 
radius of perforation are determined. Certain objections can 
be raised, however, to the results of reference 5* The 
expressions for velocity, displacement, strain, and strain 
rate presented in reference 5 <3-° not correspond to the results 
presented earlier in references 2 and 3 as expected when the 
yield strength is taken to be zero. Furthermore, the solution 
of reference 5 yields the erroneous result that for the lim- 
iting case of an infinite mass projectile the velocity of the 
projectile decreases after impact. These considerations have 
provided the motivation for the present analysis. 

The present solution agrees with the solution presented 
in references 2 and 3> in the limiting case when the yield 
strength is set equal to zero. In addition, the expression 
for velocity limits properly, for the case of the infinite 
mass projectile, to the condition that the velocity of the 
projectile remains constant after impact. The solution of 
reference 5> when compared with the present solution, is shown 
to contain additional mass ratio terms not found in the pres- 
ent solution. 

The present report compares results obtained by including 
yield strength in the analysis with the results obtained for 
the same conditions when yield strength is taken to be zero. 
Some pertinent parameters such as initial velocity, mass, and 
radius of the impacting projectile as well as the target plate 



thickness, yield strength, and viscosity axe varied in order to 
study their effects on the radius of perforation. The ballis- 
tic limits of thin target plates are investigated as well and 
the importance of yield strength considerations is ascertained. 

The ballistic limit studies determine the minimum target 
plate thicknesses necessary to prevent perforation. Ballistic 
limits are obtained for various values of target plate thick- 
ness and projectile mass, velocity, and radius. 

SYMBOLS 

The units used for the physical quantities defined in this 
paper are given both in the International System of Units (Si) 
and in the U.S. Customary Units. Factors relating the two sys- 
tems are given in reference 14 - . 
a radius of projectile 

C 1 dynamic ultimate yield strain in shear 

g Q initial projectile velocity 

h plate thickness 

i &tah/M 

H 2 rta?hp/M 

i n erfc x nth integral of erfc x (See appendix, eq. (A38)) 
I n modified Bessel function of first kind of order n 

J n Bessel function of first kind of order n 

k dynamic yield stress in shear 

K ka/V 0 |i, Bingham number 

K n modified Bessel function of second kind of order n 

Ml 


mass of projectile 



M 


r 

r 

s 

Sn 

t 

t 

V 

V 
Vo 

v_ 

Vo 

Sw _a_ 
Sr V Q 

w 

w _SL_ 
a 2 V 0 

Sw _L_ 

Sr aV Q 


mass of projectile and plug of plate material of 
radius a, Mp + jta 2 hp 
radial distance 

nondimens ional radial distance, r/a 

transform parameter 

Struve function of order n 

time 

nondimensional time, — t 

aP 

axial velocity- 
transform of V 

initial velocity of projectile and target plug of 

M-, 

radius a, V 0 = g Q -± 

M 

nondimensional velocity 

nondimensional shear strain rate 

axial displacement 
nondimensional displacement 

nondimensional shear strain 


Yn Bessel function of second kind of order n 

7 shear strain, 

Sr 

P coefficient of dynamic viscosity 

v coefficient of kinematic viscosity, jj 

p mass density of plate material 

T rz transverse shear stress 


z 


axial coordinate 


I 


Subscripts: 
cr critical 

p perforation 

A dot over a symbol denotes differentiation with respect 
to time t . 


ANALYSIS 


Governing Equations 

In the present analysis, a rigid circular cylindrical 
projectile is considered to impact upon a thin infinite plate. 
The resulting perforation of the plate by the projectile is 
assumed to be a simple shear plugging perforation, in which 
only the transverse shear stresses act to resist the inertia 
of the impacting projectile. The perforation is also con- 
sidered to be axially symmetric and the shear stress is taken 
to be constant through the thickness of the plate. The 
resulting deformations w of the plate are then represented 
as functions only of the radial coordinate r and the time t 
being independent of both the axial coordinate z and the cir- 
cumferential coordinate 0. (See fig. 1(a).) Hence, by 
summing the forces in the axial direction on a circular ring 
element of plate material (see fig. 1(b)) the basic equation 
of motion for this simple shear perforation model cam be 
written as 

^ T rz T rz 

~ + r = P 


( 1 ) 


This approach and the following discussion on plate mate- 
rial and governing equations is similar to that presented in 
references 1, 2, 3, and 6. 

The plate material considered herein is assumed to behave 
like an incompressible, visco-plastic Bingham solid in which 
the deformation accompanying transverse shearing commences only 
after the dynamic yield strength of the material has been 
reached. (See ref. r J.) When the value of the transverse 
shear stress falls below the yield strength of the material 
(or when the local rate of deformation becomes equal to zero), 
visco-plastic flow ceases and the material is assumed to be 
rigid in that region. The relation between the shear strain 
rate and shearing stress for the case of simple shear perfora- 
tion can be written as (see ref. j): 


' rz 


- •* g ♦ k i* 


(2a) 


or, since the sign dw/dr is always negative in this case 


dw 


T r* = “ £ - k 


(2b) 


where 

p dynamic viscosity of target material 

k dynamic yield strength in shear of target material 

^w Sv 

— = shear strain rate 

dr or 

and the relations between the velocity V, aocial displace- 
ment w, and shear strain 7 are 



dr 

Substitution of equation (2b) into equation (l) yields the 
governing linear partial differential equation 

d^V + ldV_ldV_lk 
3 r 2 r dr v dt r p 

where 

n 

v = ^ coefficient of kinematic viscosity 
The initial conditions are taken to be 


V = 0 

(t = 0, r > a) 

00 

V = V 0 

(t = 0, r % a) 

( 5 ) 


where V 0 is the initial velocity of the combined mass of the 
projectile and plug of target material of radius a. In this 
study V Q is determined from the conservation of momentum at 
the instant of impact in which the projectile-plug combina- 
tion is assumed to be rigid and to act as a unit. Thus 

M l 

v o - 60 M- < 6 > 

where 

g Q initial velocity of free projectile 

Mp mass of projectile 

M = Mi + «a 2 hp combined mass of projectile and plug 
The boundary conditions are as follows: at t > 0, r = a 


from Newton's second law 



where 
t s 

^ M 

and at t ^ 0, r -> 00 

v -» o (8) 

Method of Solution 

An analytical solution to the governing equations is 
obtained in reference 1 by Laplace transform techniques in a 
manner similar to the approach used in reference 5* The homo- 
geneous solution can be written directly in terms of modified 
Bessel functions. The particular solution, however, is obtained 
by utilizing Struve functions of order zero rather than by the 
method of variation of parameters used in reference 5* The 
solution to the governing differential equation, equation ( 3 ), 
and its associated boundary conditions (eqs. (4), (5), (7)> 
and (8)) is discussed in detail in the appendix. The solutions 

obtained are "short time" solutions in t corresponding to 

2 

values of the transform parameter s » ~ or t « (see 
appendix) resulting from the use of asymptotic expansions 
in the transformed state. It is shown in the appendix that 
the velocity (A3^+), displacement (A 37 ).> shear strain (A36), 
and shear strain rate (A35) can be expressed in dimensionless 
form in terms of a dimensionless radius r and a dimensionless 
time t. 

When the yield strength k is taken to be zero in the 
equations given in the appendix they reduce to the identical 
expressions given in reference 3. Also in the case of the 


* 


infinite mass projectile in which H = 0 (Mi = <») the expres- 
sion given in equation (Aj4) for velocity reduced to V = V Q 
at r = 1 for all values of t •= 0. 

Separation Criterion 

In order to determine the radius of the plug of plate 
material which is sheared (or perforated) from the plate upon 
impact a separation criterion (identical to that discussed in 
refs. 1, 2, 3, and 5) consisting of two conditions is estab- 
lished. The first condition for separation is based on the 
assumption that detachment of the projectile and a portion of 
the plate material occurs when the plate material can no 
longer transmit a given shear stress. If the plate behaves in 
a highly viscous manner it is capable of transmitting a shear 
stress of any magnitude. In this analysis it is assumed that 
the material is highly viscous when the strain rate is above 
some limiting value, here assumed to be k/p. Below this 
value only a limited shear stress can be transmitted and con- 
sequently separation may take place. This condition of high 
shear stress and visco-plastic flow will exist immediately 
after impact when the strain rate is at a max imum and the 

viscous stress p ^ is much greater than the yield strength 

dr 

k of the material (see eq.. (2b)). In equation form the 
strain-rate criterion for separation is 




The second condition deads with the shear strain of the 


material which, at the moment of impact, is zero and increases 
thereafter. In order for the plug to separate from the plate 
the material not only has to he considered as having a suffi- 
ciently small strain rate but must experience large shear 
strains as well. In this analysis the shear strains are con- 
sidered large when the shear strain exceeds the dynamic value 
of the ultimate shear strain of the material. 

Thus the second condition can be written as 



where = dynamic ultimate strain in shear (a value of 
Ci = 0.02 was used in this paper for both aluminum and steel). 
For a further discussion of these separation conditions, see 
references 1, 2, and 3* 

The radius of separation of the plug material and the 
time of separation can most readily be determined by a graphi- 
cal procedure. In figure 2 a typical nondimens ional set of 
strain distributions for different values of the time param- 
eter t is shown. In figure 3 nondimensional strain-rate dis- 
tributions are shown for a projectile of mass 0.052 g and ini- 
tial velocity 1.2 km/s impacting on a 0.635-cm (l/4-in.) alu- 
minum plate. A plot of the locus of critical strain points in 
a r,t space can now be obtained from figure 2 by constructing 
a horizontal line representing the critical strain (eq. (10)) 
and determining its r,t intercepts with the strain curves. 



In a similar manner the critical strain- rate locus in r,t 
space can be determined from figure 3* The radius of perfora- 
tion and the time of separation are now determined from the 
intersection of the critical strain and critical strain-rate 
curves as plotted in figure k. Thus for this particular case 
the radius of perforation and the time of separation as found 
from figure k are ?p = 1.7^0 and t = 0 . 250 , respectively. 

Note that four regions are indicated in figure 4. In 
these regions the predominant material behavior could be 
described as follows: large strain rates, large strain 

(region l)j large strain rates, small strain (region II); 
small strain rate, small strain (region III); small strain 
rate, large strain (region IV). Only in region IV are both 
critical conditions of the separation criterion satisfied, and 
since separation will occur at the first time at which both 
conditions are satisfied it is assumed that at the apex of 
this region separation occurs. After initial separation the 
solution is no longer valid due to the changes in the boundary 
conditions . 

EFFECTS CEF PERTINENT PARAMETERS ON PERFORATION RADIUS 

Expressions are derived in the appendix for the velocity, 
displacement, shear strain, and shear strain rate written in 
nondimensional form. These expressions contain pertinent 
parameters associated with the perforation problem such as 
projectile mass, radius, and initial velocity, and target 



density, viscosity, yield strength, and thickness. Variations 
in these pertinent parameters are now investigated to deter- 
mine their effects on the perforation radius. Numerical 
results are obtained by use of a high-speed digital computer. 
Perforation Radius as a Function of Initial Projectile 
Velocity 

In figure 5 is shown the variation in perforation radius 
versus initial projectile velocity for a 0.635-cm (l/4-in. ) 
aluminum plate and a 0.224-cm (0.088-in.) steel plate being 
impacted by a 0,119- cm (3/64-in,) radius rigid cylindrical 
projectile. The mass of the projectile was taken to be 
O.O52 g which for aluminum corresponds to a projectile 
0.434 cm (0.I7I in.) long. Note that the mass per unit area 
of the steel plate is equivalent to the mass per unit area of 
the aluminum plate. Thus the momentum exchange between the 
impacting projectile and both the steel and aluminum plates 
is identical. 

As can be seen from figure 5j the inclusion of material 
yield strength in the analysis has little effect on the perfo- 
ration radius for the range of velocities shown (maximum dif- 
ference is of the order of 5 percent). Although steel has a 
higher yield strength than aluminum, the critical strain rate, 
which depends on viscosity as well as yield strength, is for 
steel, only 2/3 of that for aluminum ^=15 ^ S - (0.15 mega- 
poise), k = 1.38 (200 ksi) for steel; p = 5 (0-05 mega- 

poise), k = O.69 (100 ksi) for Alj . Hence, for the same 


initial velocity, the perforation radius in steel is larger 
than in aluminum and separation occurs in the steel plate 
after a greater passage of time. 

Perforation Radius as a Function of Projectile Mass or 
Plate Thickness 

The effect of varying the thickness of the aluminum plate 
on the nondimensional perforation radius is shovn in fig- 
ure 6. The initial velocity of the aluminum projectile was 
taken to be 6.1 km/s and the plate thickness was varied from 
0.635 cm to 6.35 cm (l/^+ in. to 2.5 in.) all other param- 
eters being held constant. The increase in plate thickness 
decreases the perforation radius and time of separation. The 
inclusion of yield strength in the analysis has little effect 
on the calculated perforation radii (maximum variation in cal- 
culated radii is of the order of k- percent ) . Note in addition 
if the mass of the projectile is decreased by a factor (the 
radius remaining constant) instead of the plate thickness 
being increased by the same factor, the calculated results are 
identical. This result is due to the fact that the nondimen- 
sional parameter H, 


H = 


1 + 


fefr) 


(U) 


na^p 


changes in the same proportion with decrease in or 

increase in h. The perforation radii are thus seen to 
decrease with either a decrease in projectile mass or an 
increase in plate thickness. It is evident from figure 6 



that the effects of yield strength need not be considered in 
the calculation of perforation radii. 

Perforation Radius as a Function of Coefficient of 
Dynamic Viscosity 

The variation in nondimensional perforation radii due to 

changes in the coefficient of dynamic viscosity is shown in 

figure 7, for the 0.052-g projectile impacting on the aluminum 

plate. A velocity of 1.2 km/s was chosen and the coefficient 

of dynamic viscosity was varied from 0 to 15 —5— (0.15 mega- 

m^ 

poise ), all the other parameters remained constant and are 
given in figure 7* Although calculations were made for larger 
values of p, it was found that the dimensionless separation 
times exceeded the limits of the applicability of the short 
time solution. 

Note that as the coefficient of dynamic viscosity is 
increased, the resulting perforation radii also increase. A 
small variation in the coefficient of viscosity is shown to 
produce a large variation in perforation radii, and indicates 
a definite need for better evaluation of this coefficient. 

The inclusion of the yield strength in the analysis results in 
perforation radii as much as 10 percent lower than the solu- 
tion with yield strength taken to be zero. However, this 
variation can be considered to be negligible in comparison to 
the variation caused by changes in viscosity. 



Perforation Radius as a Function of Projectile Radius 


The effect of varying the projectile dimensions on the 
perforation radius (in units) is shown in figure 8. The mass 
and impacting velocity of the projectile were considered to 
remain constant at 0.052 g and 6.1 km/s, respectively. The 
target plate was considered to he made of aluminum of 0. 655- cm 
(l/ 4 -in.) thickness. The figure illustrates that as the 
shape of the cylindrical projectile changes from that of a 
pencil to that of a disk (note projectile schematics in 
fig. 8) the perforation radius increases. It is important to 
recognize that the ratio of perforation radius to projectile 
radius is diminishing with increasing projectile radius. The 
effect of including dynamic yield strength in computing perfo- 
ration radii is seen to be negligible in comparison to the 
effects produced by varying the projectile radius. 

Perforation Radius as a Function of Dynamic Yield Strength 
for a Fixed Separation Criterion 

The calculations presented earlier for the aluminum tar- 
get plate have assumed a value of dynamic yield strength in 

shear of 0.69 (100 ksi). A great deal of uncertainty 

m^ 

exists in determining the dynamic yield strength especially 
under hypervelocity Impact conditions. To determine the 
effects of wide variations in strength on the computed perfo- 
ration radii the dynamic yield strength has been varied from 0 
to 6.9 (1 000 ksi) and the calculated perforation radii are 

m2 

shown plotted in figure 9. These perforation radii are 



determined using a constant value of the critical strain rate 

= O.0692V If the critical strain rate, (k/|i), 
dr/cr *0 / 

has also been varied in accord with the variation in dynamic 
yield strength k, keeping p. constant, then as the yield 
strength approaches zero the critical strain rate would also 
approach zero and the perforation radius would approach 
infinity. This result is unrealistic and the critical strain 
rate therefore is not varied according to the variation in k 
but is assumed to remain constant. 

Figure 9 again illustrates that the effect of including 
dynamic yield strength in the calculation of perforation radii 
is negligible for reasonable values of dynamic yield strength. 
As can be seen from the figure, the decrease in perforation 
radii for an increase of the dynamic yield strength between 0 
and 5.45 ^75 (500 ksi) is less than 10 percent. 

Comments on the Effects of Yield Strength on Velocities, 
Displacements, and Stresses 

In the previous sections, it has been shown that the 
results for perforation radii obtained with the inclusion of 
yield strength in the analysis differ little from the results 
in which yield strength was taken to be zero. In reference 1 
many calculations have been made on velocities, displacements, 
and stresses with and without yield strength considerations 
and a few pertinent remarks on these results are in order. 

The strain and strain rates are little affected by the 
inclusion of yield strength terms and result in similar curves 



when plotted as a function of plate radii. The velocities, 
displacements, and stresses of the solution containing yield 
strength, however, are at times markedly different from the 
solution in which yield strength is taken to be zero. The 
effects of yield strength in general are to produce much 
smaller velocities and much smaller displacements than when 
the yield strength is taken to be zero. In fact, for the par- 
ticular case where the projectiles do not perforate the target 
plate, the analysis including yield strength indicates the 
velocity drops to zero and the plate becomes rigid in a 
deformed state, while if the yield strength is taken to be 
zero the velocity approaches a nonzero limiting value, and the 
plate continues to deform. A more complete discussion and 
analysis of the effects of target material yield strength on 
the velocities, displacements, and stresses of the target plate 
is given in reference 1. 

BALLISTIC LIMITS OF THIN PLATES 

Another question of prime concern in the investigation of 
hypervelocity impact is that of establishing mini mum thickness 
requirements of thin plates necessary for the prevention of 
projectile perforation. This minimum plate thickness or maxi- 
mum projectile velocity allowable before perforation initiates 
is termed the ballistic limit. The ballistic limit was deter- 
mined in the following manner. The graphical procedure as 
illustrated in figure k by the intersection of the critical 



strain and critical strain-rate curves, establishes the radius 
and time of separation of the plug of plate material being 
sheared from the target plate. By decreasing the Initial pro- 
jectile velocity, the critical strain curve in figure 1+ rises, 
and the critical strain-rate curve decreases and bends more 
sharply toward the ordinate t. Thus, the radius of perfora- 
tion and time of separation decrease until the intersection of 
the two curves occurs at r = 1. For any further decrease in 
initial projectile velocity, the curves will not intersect and 
a ballistic limit is assumed to have been reached. Note that 
under the assumptions adopted in this paper no considerations 
have been made for spallation. The ballistic limits estab- 
lished in this manner for thin target plates could yield plate 
thicknesses which are unconservative (i.e., too thin to pre- 
vent perforation for a given projectile velocity) and compari- 
son with experimental data is essential to ascertain the 
degree of error introduced in neglecting the effect of spalla- 
tion of the target plate. 

The ballistic limit thickness as a function of projectile 
mass for six different projectile velocities is shown in fig- 
ure 10. The projectile velocities were varied from 0.91^ km/s 
to 6.462 km/s. However, the projectile radius was held 
constant at 0.32 cm (l/8 in.). The curves which are shown in 
the figure are valid for k = O.69 — (100 ksi) and the projec- 

m2 

tile velocity shown above the curves, and for k = 0 and the 
projectile velocity shown below the curves. The inclusion of 



the material yield strength in ballistic limit studies produces 
effects which are only slightly higher (approximately 10 per- 
cent) than those presented earlier on perforation studies 
(approximately 5 percent). 

The minimum plate thicknesses are shown plotted as a 
function of initial projectile velocity for three different 
mass projectiles in figure 11. Figure 11 is simply a cross 
plot of figure 10, considering the mass Mp constant. The 
radius of the projectile is held constant at 0.32 cm 
(l/8 in.), but the mass is varied from 0.001 g to 0.100 g. 

This change in mass could be accomplished by either changing 
the length or the density of the projectile. Note that the 
curves are all parallel and increase as (g 0 )^*®^. This 
velocity dependence is quite similar to experimental results 
obtained for the- dependence of penetration on the velocity of 
the impacting projectile by Collins and Kinard (ref. 8). 

In figure 12 the minimum plate thickness is again plotted 
as a function of projectile mass for eight different projec- 
tile velocities. The radius of the projectile, however, has 
now been decreased to 0.119 cm (3/64 in.). Notice how 
sharply the minimum plate thicknesses have increased from fig- 
ure 10 to figure 12. In the case of a projectile of mass 
0.01 g, the minimum plate thickness differs by a factor of 
five and a half considering the projectile first to have a 
radius of 0.32 cm (l/8 in. ) and then a radius of 0.119 cm 
(3/64 in.). Under these conditions, the same density 



material would resemble first a disklike projectile and 
then a bulletlike projectile as the radius is decreased. 
Hence, projectile aspect ratio is most important in deter- 
mining the ballistic limit of thin plates, exhibiting a much 
more influential effect on the results than is evident in 
perforation studies. 


CONCLUDING REMARKS 

It is shown that the inclusion of target material yield 
strength in a one-dimensional analysis of hypervelocity impact 
perforations produces little effect on the resulting perfora- 
tion radius. The variation in calculated perforation radius 
as compared to the solution with target yield strength taken 
to be zero generally amounted to only 5 percent in this inves- 
tigation. This difference in calculated perforation radii 
remains at the negligible 5-percent level even with variations 
in the pertinent parameters such as plate (or target) thick- 
ness, projectile mass, and initial velocity. With variations 
in the coefficient of dynamic viscosity, however, the differ- 
ence due to yield strength does increase to perhaps 10 percent. 
On the other hand, the variation in perforation radius due to 
differences in the assumed value of the viscosity coefficient 
alone can be much greater than the 10 percent. Consequently, 
the determination of accurate values for the coefficient of 
dynamic viscosity is much more critical to the calculation 
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of perforation radius than is the inclusion of the target 
yield strength. 

In the determination of the ballistic limits of thin 
plates (that minim um allowable thickness or maximum allowable 
velocity) the effect of including yield strength in the analy- 
sis produces results which differ by approximately 10 percent 
from the analysis with yield strength taken to be zero. A 
much more critical parameter was found to be the radius of 
the impacting projectile. The aspect ratio of the Impacting 
projectile (in this case a right circular cylinder), in 
changing from a disklike projectile to a pencillike projectile 
of the same mass, produces ballistic limits which differ con- 
siderably. Obviously the radius of the impacting projectile 
is therefore a factor of prime importance, and the yield 
strength of the plate material one of secondary importance in 
ballistic limit studies. It was noted that the nearly linear 
velocity dependence in the determination of ballistic limits 
was quite similar to experimental results obtained in the 
penetration studies of reference 8. 



APPENDIX 


GOVERNING EQUATION OF SHEAR PERFORATION 


In this appendix an analytical solution is presented for 
the governing linear differential equation and its associated 
boundary conditions. In a manner similar to that employed 
by Chou in reference J>, the governing partial differential 
equation is reduced to a total differential equation by use of 
Laplace transform techniques. To obtain a particular solution 
to the resulting total differential equation a further change 
in variables becomes necessary. After satisfying the asso- 
ciated boundary conditions a "short time" solution is deter- 
mined by employing asymptotic approximations in the trans- 
formed state. 

The governing linear differential equation of motion in 
the axial direction, for simple shear perforation, as derived 
in the text (see eq. ( 3 )) can be written as 


d 2 V + ldV_ldV = lk 
d r 2 r dr " v 5t r |i 

transforming equation (Al) with respect to t by use of 
Laplace transform techniques, results in 


(Al) 


— + 
dr^ 



1 k_ 
r |is 


(A2) 


where s is the transform parameter. By use of the initial 
condition (eq. (l^)), V(o,r) =0 at t » 0, 
equation (A2) is reduced to 


r > a 



(A3) 


afv + ldV_sV_lk_ 
dr2 r dr “ v ' r ns 


Let 


-iz 



(A) 


hence 


dV _ dV . £ 

dr “ dz VV 

d^V _ -d 2 ? ^ 
dr 2 dz 2 v 

and by substitution of the above relations into equation 
there results 


(A5) 

(A6) 

(A3) 


d 2 ^ , ldV^.^LJ. 
dZ 2 Z dZ yjs ^ sz 


(AT) 


The homogeneous solution can he written directly as 

V H (z) = AJ 0 (z) + BY 0 (z) (A3) 

For the particular solution to equation (A7) the Struve func- 
tion of zero order is satisfactory. A known relation for the 
Struve function of order p is (see ref. 9, p. 211) 


d 2 Sp(z) i dSp(z) 
dz 2 + z dz 




\/jc(p - 0.5)! 


(A9) 


Let p = 0 


d 2 S 0 (z) -l dS Q (z) 

+ 

dz 2 z dz 



(A10) 


Hence, the particular solution to equation (Aj) can be written 


as 


V p (z) 


2Si_ 


ff i 


ns c 


(z) 


(All) 


and the general solution becomes 


V(z) = V H (z) + V p (z) = AJ 0 (z) + BY 0 (z) 2^_S 0 (z) 

2 I¥ S 


or, resubstituting for z from equation (Ah) 

V(r,s) - iJ„(ir {£) ♦ BY 0 ( lr /I) - 0 (lrji) 


(A12) 


(A13) 


In order to apply now the (transformed) boundary condition 
given by equation (8), namely, V = 0 as r -* oo one must con- 
sider the asymptotic behavior of the solution as given by 
equation (A13) for |r| » 1. The asymptotic approximation for 
the Struve function of order zero as found in reference 9, 
equation (136a) is 


®o 




(AIM 


and since a known identity between the Bessel functions is 
(see ref. 10, eq. (110)) 


*°( r if) - l 1 ^ if) - 1Y °( lr if) 

then 


(A15) 



(A16) 


Hence, for r | >> 1, the general solution can be written as 



V(r,s) - AI 0 (r jf ) + B il 0 (r JI) - )j 


+ 




(A17) 


Application of the transformed boundary condition, V = 0 as 
r -»« (see eq. ( 8 )) yields the relation 


A « -B ± - 


jtk 


2 (A18) 

Since Kq -* 0 and I Q » as r -»oo. Use of the known iden- 
tities between Kq (rv/|) (eq. (A15)), S G (ir (eq. (Al4)), 

X 0 (r$, and Y 0 ( ir 4 ) reduces the general solution to 
V (r ,s) , -B^(r £) - -f-! 0 (r jf) - -^o(ir ft) (A19 , 






The quantity B is now evaluated by application of the 
(transformed) boundary condition (see eq. (j); note also 
eq. ( 5 )): 


+ C| - V 0 (r - a, t > 0) (A2Q) 

Substitution of equation (A19) into equation (A 20 ) and by use 
of equation (l40a) of reference 10 yields the following expres- 
sion for B: 





Hence, the general solution can he written as 



A "short time" solution is now determined from the general 
solution by assuming a ^ > 1. By reference again to equa- 
tion (lj6a) of reference 10, it can be shown that 
S.l(ia la #) + | 2 




(A23) 


A known identity between the modified Bessel function of the 
second kind of order one K]_ and the Hankel function of the 
first kind of order one is (see ref. 10, eqs. (115) 

and (70)) 


% ( a /I ) = - i H i 1) ( ia /!) - - § (ji (** jl) + (*» /?) 


(A2k) 


Hence 


j i( ia Jl)- 

and since 

x i( a Jv) = _iJ i( ia /f) 


(A25) 


(A26) 


then, from equation (A25) 


Il(a jf) » /|) + §*l( a /|) (A27) 


t - s l( la if) ~ J l( a /I) + - f*l( a Jf) 


(A28) 


Substitution of the asymptotic expressions for S Q and 
^ - Sq into the general solution (eq. (A22)) and making use 
of the identity for Y 0 (eq. (A15)) yields 


V(r,s) = 


( v ° + £s - S ♦ ^) K °( r /?) 

■*»(» if) + c* iK( a if) 


(A29) 


Note that for increased accuracy in equation (A29), it is 
only necessary to increase the number of terms taken in the 
asymptotic expansions of S 0 and S^; of course, the semicon- 
vergent nature of the asymptotic expansions must be considered. 
These additional terms are polynomials and would not present 
any new difficulties. To aid in determining the inverse trans- 
form of equation (A29), the following asymptotic expansions of 
the modified Bessel functions Kq, K]_ are employed (see 
eq. (114) of ref. 10) 


M*) - - £ 


+ • • • 

8z 128z 2 


K l< z > - + £ 


128z 2 


(A30) 


Substitution of the asymptotic expansions given in equa- 


tion (A30) into equation (A 2 $) yields, after algebraic 
manipulation 


V(r,s) - ^ - 5t) * ^][l - ($ * 

+ ( 2 *- Z 2ar ~ 111 + ft a - 5 r f 2 ^ 1 "1 k 

V 128 a 2 r 2 J v 8ra 5 v 2 /^ J “ ^ (AJl) 


9a £: - 2ar - 7r^ , ^ a - 5r 


+ ^ is + tiS s, - 2* + 

V 128a 2 r 2 v 8ra 


and by application of standard transform tables (ref. 11, 
p. 380, formulas (ll)) V becomes 

V(r,t) = V 0 (a) l/2 Lfc JL^_a + (£-=_a - erfc r^ji 

W L 2j^ V 8r W a sfvt 


-- 7r2 ♦ # ■ ■ ] 

(= - (Hr - ^ 

9a 2 - 2 ar - 7r 2 pa a - 5r 2 p 2 a 2 \l6vt 2 ., r - a 1 

128,2 ? v a, 5 -?r)-sr - 1 ' ' J 


9a 2 - 2ar - 7r g 


128 r 2 


^(|W 2 ri6t 2 i^erfc + (LJl* - ^ 

pa 2 V r / [_ 2 f t V 8r 


9a 2 - 2ar - 7r 2 
1 128 r 2 


+ ^ a-5r + ^2 H^WtVWfc 


v 8r 


V 2 /a 2 


i kt 
J ” pr 


where the symbol i n denotes the nth integral of the comple- 
mentary error function. 

Equation (A32) can now be nondimens ionali zed as follows: 


Let 


— r 
r - a 


f _ vt 


t = 

2 > 

n 2«a lip 

^ V M 


K = 


ka 

V 


(A33) 


Hence, the nondimens i onal velocity can be written as 


i = 4 <srtc r * 1 






+ 1«(X - H)tl 2 erfc - ~_} + XfiKHi^erfc r " 1 


2<ft 


2^ 


erfc - -_J- + GK(l - H)tl5 e rfc r ' 1 + 32KH(t) 5 / 2 l 5 erfc 

2 VT 


2$ 


V] 


+ RxJWti 

+ R 2 l'+tl 2 erfc --- 1 - + 16K(1 - H^lLrfc r -~ ^ + 64KfiF 5 erfc r -~ - 1 ~ |. . i - £[ 

L 2& ^ 2\l7 J J r 

((F - l) £ o) 


(A 34 ) 


where 



JL 

8r 


H 


R 2 


128F 2 


JL 7 

6kr m 




+ H 2 


By comparison of equations (A 3 l) and (A32), it can be shown 
that s is related to l/ 4 t and that a^s/v > 1 implies 
t < 1 / 4 . Nevertheless, calculations indicate that the three- 
term asymptotic series presented in equation (A30) are only 
approximately 10 to 15 percent in error for t = l/2 (and 


\ 



a^/V = 1/^2 j. Therefore, results are presented in this paper 
for values of t as high as l/ 2 . Results for t > l /2 should 
he viewed with increasing skepticism. Differentiation of equa- 
tion ( A38) with respect to r yields the nondimensional shear 
strain rate 


(r-l)2 


+ R3 erf c + 2jt(R4 - K(l - H)] i erfc 


Sr V ° ^F Jrtt J 2 ft L ^ 2 \|t 

+ it[R 5 + K(l - H)R^ji 2 erfc r ~ 1 + 8(t) 5 / 2 [k(l - H)Rl* - KH]i 5 erfc r ~_3 

2ft 2ft 

+ l6t2[KHR 5 + K(l - H)r£J i^erfc £ yl + 32(t) 5 / 2 KHR4i 5 erfc r 


+ eit 5 ™^ 6 erfc 1,-1 . . A + % 
5 ( ? 


((? - 1) 5 0) 


where 


(A35) 


o _ 3 1 . - 

5 "8r " 8 H 


R 4 = 


_ 15 


5 + J- + 


I28F 2 6 k¥ 128 


If 5 )- 


R c 


Rr = - 


^ + 1 + JL_ + ML 

256r 5 128T 2 2 56r 16 \r 


£ 

2r 


By integration of equation (A35) with respect to t the non- 


dimensional. shear strain is obtained 


2H. _v_ = erfc r ~3 + R^i 2 erfc r 1 + 2^fR4 - K(l - H)1 inerfc r - 1 

dr V o a iF][ 2 Jf 2 ^ 5 2 /t L J 2 jt 

+ 4t[R 5 + K(l - H)r£J inerfc f ~ -i + SCt)?/ 2 ^! _ H )R 4 - KH]i 5 erfc ? -J 1 - 

2/f 2/t 

+ l^fKHRj + K(l - H)R£Ji 6 erfc f ~- A + 32(t) 5 / £ KHRi l .i’ 7 erfc £ - ~ 1 

2ft 2W 


+ 6Ut 5 KHR 5 i 8 erfc r "A . . A + — 


2 it" f S ? 2 


((7 - 1 ) £ 0) 


(A36) 


A second integration with respect to r would determine the 
corresponding displacement. However, a simpler approach is to 
integrate the nondimensional velocity with respect to t to 
obtain the nondimensional displacement 


w-^— = it! /i 2 erfc + 4K(l - H)tl**erfc £- =- 1 + 16 KHt^erfc r ~ 1 

a V ° 2||7 2 U 2 \T 

+ R x J2 if i ^erf c ? + GK(l - H) ( t ) 5 / 2 i 5 erf c E-=_i + 32KH(t) 5//2 i 7 erfc 

L 2l|f 2 if 2 if J 

[utiAerfc r + l6K(l - Hjt^erfc r ~_ X + 64KHpi®erfc . . \ - 

L 2 if 2jt 2jt J j 


+ R 2 

L 2l|t 2jt 

Vr - 1)i0> (A57) 

The symbol i n erfc x represents the repeated nth integrals of 
the complementary error function erfc ! 


i n erfc x 


J ^OO 


erfc| d! 


(A38) 


,0 


erfc x - erfc x 


i -1 erfc x = ie- xC 

J” 


in which 



A recurrence relation useful in determining the repeated inte- 
grals of the complementary error function is (see refs. 10 
and 11) 


i n erfc x - 


i n ~ 2 erfc x - 2xi n ~ 1 erfc 


A table of the error function and its derivatives and integrals 
for values of the argument between 0 and 3*0 is presented in 
reference 11 from values given in reference 12. A similar 
table has been included in reference 1 with slight extensions, 
for convenience of use in computing the above expressions. A 
table of the first 11 repeated integrals of the error 
function with values of the argument between 0 and 2.80 is also 


available in reference 13 . 
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Figure 1.- Geometry and coordinate system of perforation. 
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Figure 5*- Radius of perforation as a function of initial projectile velocity for steel and 
aluminum target plates of the same weight per unit area. Mp = 0.052 g; a = 0. 1 I Q cm 
(3/64 in.); ph = 17 £§•; p = 5 —• (0.05 megapoise). 
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, VARIES WITH COEFFICIENT OF DYNAMIC VISCOSITY (4) 



Figure 7.- Effect of variation in coefficient of dynamic viscosity on perforation radius 

h = 0.635 cm (1/4 in.); M-, = O.O52 g; g Q = 1.2 km/sec; p = 2680 £§; a * 0.119 cm 

"" va? 

(3/64 in.). 






CRITICAL STRAIN-RATE ■ -0.069 
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Figure 11.- Minimum plate thickness as a function of projectile velocity for projectiles of three 
different masses, p = 5 (0.05 megapoise); p = 2680 £f; k = O.69 ~ (100 ksi) a = 0.32 cm 

(l/8 in. ) . 
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Figure 12.- Ballistic limit curves for projectile of 0.119 cm (3/64 in.) 
radius, p = 5 (0.05 megapoise); p = 2680 


NASA -Langley, 1964 



